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Abstract 

(— I , The completely bounded trace and spectral norms, for finite-dimensional spaces, are known 

p I to be efficiently expressible by semidefinite programs (J. Watrous, Theory of Computing 5: 11, 

2009). This paper presents two new, and arguably much simpler, semidefinite programming 
formulations of these norms. 

1 Introduction and preliminary discussion 

In the theory of quantum information, quantum states are represented by density operators act- 
ing on finite-dimensional complex vector spaces, while quantum channels are represented by linear 
mappings that transform one density operator into another I NC00llKSV02 l. Various concepts con- 
nected with mappings of this form, meaning ones that map linear operators to linear operators (or, 
equivalently, that map matrices to matrices), are important in the study of quantum information 
for this and other reasons. Linear mappings of this form are also important in the study of operator 
algebras [Pau02J. 

This paper is concerned specifically with the completely hounded trace and spectral norms, defined 
for linear mappings of the form just described. It is intended as a follow-up paper to [Wat09'l, 
which demonstrated that these norms can be efficiently expressed and computed through the use 
of semidefinite programming. Two new semidefinite programming formulations of these norms 
will be presented, both of which are simpler than the formulations given in the previous paper. 

A further discussion of the completely bounded trace and spectral norms can be found in 
|jWat09]. That discussion will not be repeated here — instead, we will proceed directly to the tech- 
nical content of the paper, beginning with a short summary of the notation and basic concepts that 
are to be assumed. 



Linear algebra basics 

For a complex vector space of the form X = C" and vectors u = [ui, . . . ,u„) and v = {vi, ...,??„) 
in X, we define the inner product 

n 

{U,V) = ^M-Z7y 
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as well as the Euclidean norm 




For each ; G {1, . . . , n}, the vector Cj G X is defined to be the vector having a 1 in entry and for 
all other entries. 

Given two complex vector spaces X = C" and y = C", we denote the space of all linear 
mappings (or operators) of the form A : X ^ y as L{X,y), and identify this space with the 
collection of all m x n complex matrices in the usual way. For each pair of indices we write 
E, y to denote the operator whose matrix representation has a 1 in entry (f, /) and zeroes in all other 
entries. The notation L{X) is shorthand for L{X,X), and the identity operator on X, which is an 
element of L(A'), is denoted Ix- (The notation 1 is sometimes used in place of lL;f when it is clear 
that we are referring to the identity operator on X.) 

For each operator A G L(X,y), one defines A* G L{y,X) to be the unique operator satisfying 
{v,Au) = {A*v,u) for aHu E X and v E y. As a matrix. A* is obtained by taking the conjugate 
transpose of the matrix associated with A. An inner product on L{X, y) is defined as 

{A,B) =Ti{A*B) 

for all A, B G L{X,y). By identifying a given vector u E X with the linear mapping a i-)- au, 
which is an element of L(C, X), the mapping u* G L{X,C) is defined. More explicitly, u* is the 
mapping that satisfies u*v = {u,v) for all v E X. 

An operator X G L(A') is Hermitian if X = X*, and the set of such operators is denoted 
Herm(A'). An operator X G L{X) is positive semidefinite if it is Hermitian and all of its eigenvalues 
are nonnegative. The set of such operators is denoted Fos{X). The notation X > also indicates 
that X is positive semidefinite, and more generally the notations X < Y and Y > X indicate that 
Y — X > for Hermitian operators X and Y. An operator X G L(A') is positive definite if it is both 
positive semidefinite and invertible. Equivalently, X is positive definite if it is Hermitian and all of 
its eigenvalues are positive. The set of such operators is denoted Pd{X). The notation X > also 
indicates that X is positive definite, and the notations X < Y and Y > X indicate that Y — X > 
for Hermitian operators X and Y. An operator p G L(A') is a density operator if it is both positive 
semidefinite and has trace equal to 1, and the set of such operators is denoted D{X). Finally, an 
operator U E L(A') is unitary iiU*\I = \x ^ and the set of such operators is denoted \J{X). 

For X = C" and y = C", the space of all linear mappings of the form $ : L(A') — > L{y) is 
denoted T(X,y). For each <I> G T{X,y), the mapping <I>* G T{y,X) is the unique mapping for 
which the equation 

(Y,0(X)) = (<J>*(Y),X) 

holds for all X G L(A') and Y G L(3^). A mapping <I> G T(A', y) is Hermiticity preserving if it holds 
that ^(X) G Herm(3;) for all choices of X G Herm(A'), positive if it holds that <[)(X) G Pos(3^) for 
all X G Pos{X), and completely positive if O (g) 1l(c*) positive for all k >1. 

Norms and fidelity 

For X = C", y = C", and any operator A E L{X, y), one defines the trace norm, Frobenius norm, 
and spectral norm as 



A 111 = Tr VA*A, \\A\\2 = J {A,A) , and || A||oo = max{|| Am|| : m G A", ||u|| < l}. 
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respectively. These norms correspond precisely to the 1-norm, 2-norm, and oo-norm of the vector 
of singular values of A. All three of these norms are unitarily invariant, meaning that 

||iiAy||i = ||LZAy||2 = ||A||2, and \\UAV\\^ = \\A\\^ 

for every operator A G L(A', 3^) and every choice of imitary operators U G U(3^) and V G \J{X). 
For every operator A G L( A", 3^) it holds that 

ll^lloo < ||A|j2 < ||A||i. 

The trace and spectral norms are dual to one another, meaning 

||A||i =max{|(B,A)| : ||B||eo < l}, 
||A||c„ = max{|(B,A)| : ||B||i < l}, 

for all A G L( A", 3^), and with B ranging over operators within the same space, while the Frobenius 
norm is self-dual. 

For each O G T( A", y), one defines the induced trace and spectral norms as 

llOlli = max{||0(X)||i : X G ||X||i < l}, 

||0||eo =max{||0(X)||eo : X G L(A'), ||X||eo < l}, 

as well as completely hounded variants of these norms: 

|||0|||i =SUp||0(g)lL(c*)||i = 
k>\ 

|||0|||oo =sup||0^1l(c*)|L = ||^^1l(3^)IL- 
By the duality of the trace and spectral norms, it holds that 

lll^llll = lll^llloc (1) 

for every mapping O G T{X,y). In the subsequent sections of the paper, our focus will be on 
semidefinite programming formulations of the completely bounded trace norm ||| ■ interested 
readers may directly adapt these formulations to ones for the complete bounded spectral norm by 
means of the relationship ([l]). As every operator X having trace norm bounded by 1 can be written 
as a convex combination of rank 1 operators taking the form uv* for u and v being unit vectors, it 
follows from the convexity of norms that 

|||0|||i =max{||(0(g)lL(;t.))(uy*)||^ : u,v e X ® X, \\u\\ = \\v\\ = l}. (2) 

Finally, for any two positive semidefinite operators P, Q G Pos(A'), one defines the fidelity 
between P and Q as 

F(P,Q)= VP^/Q . (3) 

1 

For u,v ^ X ® y being any choice of vectors, it holds that 

F{Try{mn,Try{vv*)) = ||Tr;,(i;w*) (4) 

(It should be noted that the partial traces on the left-hand-side of the equality in this theorem are 
taken over the space y, while the partial trace on the right-hand-side is taken over X.) A proof of 
this identity may be found in I.RW05 I or IIWat08H . 
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Semidefinite programming 

A semidefinite progran^ is specified by a triple (E, C, D), where 

1. H G y) is a Hermiticity-preserving linear map, and 

2. C G Herm(A') and D G Herm(3^) are Hermitian operators, 

for X = C" and y = C"' denoting spaces as before. We associate with the triple (E, C, D) two 
optimization problems, called the primal and dua/ problems, as follows: 

Primal problem Dual problem 

maximize: (C, X) minimize: {D,Y) 

subject to: S(X) = D, subject to: E*{Y)>C, 

X e Pos(A'). y e Herm(3^). 

An operator X G Pos(A') satisfying H(X) = D is said to be primal feasible, and an operator 
y G Herm(3^) satisfying E*(y) > C is said to be dual feasible. We let V and "D denote the sets of 
primal and dual feasible operators , respectively: 

7^ = {X G Pos(A') : S(X) = D} and V = {Y e Herm{y) : E* [Y) > C} . 

The linear functions X i— ?► (C, X) and Y i— )■ {D,Y) are referred to as the primal and dual objective 
functions, which take real number values all choices of X G 7^ and Y ^ D (or, more generally, over 
all choices of X G Herm(A') and Y G Herm(3^)). The primal optimum and dual optimum are defined 
as 

a = sup(C,X) and ,6 = inf {D,Y), 

respectively. The values a and (5 may be finite or infinite, and by convention we define a = —coif 
V = and /5 = 00 if D = 0. If an operator X eV satisfies (C, X) = a we say that X is an optimal 
primal solution, or that X achieves the primal optimum. Likewise, if y G I? satisfies {D,Y) = fi we 
say that Y is an optimal dual solution, or that Y achieves the dual optimal. 

For every semidefinite program it holds that a < j6, which is a fact known as zveak duality. 
The condition a = j6, known as strong duality, may fail to hold for some semidefinite programs — 
but, for a wide range of semidefinite programs that arise in practice, strong duality does hold. 
The following theorem provides a condition (in both a primal and dual form) that implies strong 
duality. 

Theorem 1 (Slater's theorem for semidefinite programs). The following implications hold for every 
semidefinite program (H, C, D). 

1. lfV^0 and there exists a Hermitian operator Y for which S* (y) > C, then a = fi and there exists a 
primal feasible operator X G V for which (C, X) = a. 

2. lfV^0 and there exists a positive definite operator X > Q for which E(X) = D, then a = (i and 
there exists a dual feasible operator Y G Vfor which (D, y) = /6. 

The condition that some operator X > satisfies E(X) = D is called strict primal feasibility, while 
the condition that some operator Y G Herm(3^) satisfies E*(y) > C is called strict dual feasibility; 
in both cases, the "strictness" concerns the positive semidefinite ordering. 



^ It should be noted that the above definition differs slightly from the one in fWat09l, where the equality con- 
straint H(X) = D appears instead as a inequality constraint E(X) < D, and (correspondingly) the dual condition 
Y 6 Herm(3^) appears as Y 6 Pos(3^) . The two forms can easily be converted back and forth, but the one above is more 
convenient for the purposes of this paper. 
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2 A semidefinite program for the maximum output fidelity 



The first semidefinite programming formulation of the completely bounded trace norm to be pre- 
sented is based on a characterization of the completely bounded trace norm in terms of the fidelity 
function, together with a simple semidefinite program for the fidelity function itself. 

2.1 A semidefinite program for tfie fidelity function 

We will begin by presenting a semidefinite programming characterization of the fidelity F(P, Q) 
between two positive semidefinite operators P, Q G Pos(A'), for X = C". The same semidefinite 
programming characterization of the fidelity was independently discovered by Nathan Killoran 

The semidefinite program is given by the triple (H, C, D), where S : L{X © A") — L{X © X) 
is defined as 



for all Xi 1, Xi^2/ X24, X2,2 S and C, D G Herm(;f ® M) are defined as 

The primal and dual problems associated with this semidefinite program may, after some simpli- 
fications, be expressed as follows: 

Primal problem Dual problem 

11 11 
maximize: - Tr(X) + - Tr(X*) minimize: -{P,Y) + -{Q,Z) 

subject to: (^^ q) " ^ subject to: (^_^ Z^) ^ ^ 



xeL(A'). y,z e Herm(;f). 



Strong duality 



Strong duality for the semidefinite program (S, C, D) may be verified through an application of 
Slater's theorem, using the fact that the primal problem is feasible and the dual problem is strictly 
feasible. In particular, the operator 

is primal feasible, which implies that V ^ 0. For the dual problem, the operator 

1 

is strictly feasible, as 

'1 0\ /I 0\ 1/0 1 



11 VO 11 ^ 2\1 



By Slater's theorem, we have strong duality, and moreover the primal optimum is achieved by 
some choice of a primal feasible operator. 
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It so happens that strict primal feasibiHty may fail to hold: if either of P or Q is not positive 
definite, it cannot hold that 

P X' 



One cannot conclude from this fact that the optimal dual value will not be achieved — but indeed 
this is the case for some choices of P and Q. If P and Q are positive definite, however, then strict 
primal feasibility does hold, and the existence of an optimal dual solution follows from Slater's 
theorem. 



Optimal value 

One may prove that the optimal value of the semidefinite program described above is equal to 
F(P, Q) by making use of the following fact (stated as Theorem IX.5.9 in IIBha97B V 

Lemma 2. Let P,Q E Pos(C") be positive semidefinite operators and let X G L(C") be any operator. It 
holds that 

^) GPos(C"©C") (5) 



if and only ifX = ^/PK^/Qfor K G L(C") satisfying \\K\ 



< 1. 



It follows from this lemma that for feasible solutions to the primal problem, the variable X G L( A:") 
(in the simplified form of the primal problem) is free to range precisely over those operators given 
by \/~PK\^ for K G L(A') satisfying ||i<^||oo < 1- The primal optimum is therefore given by 

supQTr(^KyQ) +iTr(yQK*v^)^ = sup ^(Tr (^QK* v^' 



sup 

K 



Trfygr^/p' 



sup 

K 



k,VpVq) = y/Py^ =F(P,Q), 



where each supremum is over the set {K E L{X) : ||iC||oo<l}- 

By strong duality, the dual optimum is also equal to F(P, Q). An alternate way to prove this 
fact begins with the observation that the dual optimum is equal to 



;Pd(A')V2^ '2^ 'J 



yePd(A') ' 

This expression follows from the observation that, for every Y,Z E Herm( A"), it holds that 



(6) 



y 
-1 



-1 



if and only if Y, Z G Pd( A") and Z > Y^^, together with the assumption that Q is positive semidef- 
inite. Now, the fact that the dual optimum is equal to F(P, Q) follows from a theorem known as 
Alberti's theorem. 

Theorem 3 (Alberti). Let A" = C" and let P, Q G Pos(A') be positive semidefinite operators. It holds 
that 
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To see that Alberti's theorem implies that the expression ^ is equal to F(P, Q), note first that the 
arithmetic-geometric mean inequality implies that 

for every Y G Vd{X), with equality if and only if (P, Y) = (Q, Y^^) . It follows that 

Moreover, for an arbitrary choice of Y G Pd(A'), one may choose A > so that 

(P,AY> = (Q, (AY)-i> 

and therefore 

1(P, AY) + \{Q, (AY)-i> = ^(P,AY)(Q,(AY)-i) = ^(P, Y)(Q, Y-i). 

Thus, 

inf (hp,Y) + \{Q,Y-^)]=V{P,Q). 

By reversing this argument, an alternate proof of Alberti's theorem based on semidefinite pro- 
gramming duality is obtained. A similar observation was made in IIWat09| based on a different 
semidefinite programming formulation of the fidelity. 

2.2 Maximum output fidelity characterization of the completely bounded trace norm 

Next, we recall a known characterization of the completely bounded trace norm in terms of the 
fidelity function, which makes use of the following definition. 

Definition 4. Let ^ = C" and Z = C^, and let Yq, Yi G T{X, Z) be positive maps. The maximum 
out-put fidelity between Yq and Yi is defined as 

Fmax(Yo,Yi) =max{F(Yo(po),Yi(pi)) : po,pi G D(A')}. 

The characterization (which appears as an exercise in IIKSV021 and is a corollary of a slightly more 
general result proved in [WatOSJ ) is given by the following theorem. 

Theorem 5. Let M = C", y = C", and Z = let Aq, Ai G L(A',3^ (g) Z) be operators, and let 
Yq, Yi G T{X , Z) and O G T{X , y) he mappings defined by the equations 

Wo{X) =Jry{AoXA*o), Yi{X) =Jry{A^XAl), and ^{X) =Trz{AoXAl), 

for all Xe L(A'). It holds that |||0|||i = Fmax(Yo,Yi). 

Proof. For W = C" and any choice of vectors uq, Mi G A" (8) W, one has 

Tr3;0w((^o ^ lw)woMo(^o ^ Iw)*) = Yo(Trn)(MoWo))' 
Try^wiiM (S) lw)uiul{Ai (g) IwT) = Yi {Trwiu^ul)), 
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and therefore, by it holds that 

\\Ttz{{Ao ® lw)uoul{A, ® IwT) 111 = F(Yo(Trw(wo"5))''fi(Tr>v("iwD))- 
Consequently 

IllOjlli =max{||Tr2((Ao(8)lw)wowU^i'^lw)*)||i : "O/"! e X<^yV, \\uo\\ = \\ui\\ = 1} 
= max{V {Y o{Trw {uqUq)), Y i{Jrw (uiul))) : uq,Ui E X ®W, \\uq\\ = ||mi|| =1} 
= max{F(Yo(|Do),Yi(pi)) : po,Pi S D(A:')} 

= Fmax(Yo,Yi) 

as required. □ 



2.3 A semidefinite program for the maximum output fidelity 

Theorem |5j when combined with the semidefinite program for the fidelity discussed at the be- 
ginning of the present section, leads to a semidefinite program for the completely bounded trace 
norm, as is now described. 

Let X = C" and y = C"\ and suppose that a mapping O G T(A', 3^) is given as 

cS>{X)=Trz{AoXAl) (7) 

for all X G where Z = and Aq, Ai G L(A', y (Si Z) are operators. An expression of this 

form is sometimes known as a Stinespring representation of O, and such a representation always 
exists (provided that k is sufficiently large; k must be at least mn in the worst case). 
Now, define completely positive mappings Yq, Yi G T( A", Z) as 

Yo(X) = Try{AoXA*o) and Yi(X) = TTy{A^XAl) 

for all X G L(A'). The semidefinite program to be considered is specified by the triple (S, C, D), 
where E : L{A! (B X (S) Z (B Z) — >L(C©C©Z©Z)isa Hermiticity-preserving mapping defined 
as 

■ • • \ /Tr(Xo) \ 
Xi ■ ■ _ Tr(Xi) 

■ Zo ■ ~ Zo-Yo(Xo) 
• ■ Zj \ Zi-Yi(Xi)y 

where dots represent operators on appropriately chosen spaces upon which S does not depend, 
and C G Herm(;f © A" © Z © Z) and D G Herm(C © C © Z © Z) are defined as 



V 



(8) 



= 4 



/O 0\ 



1 

Vo 1 0/ 



and 



D 



/I 













1 




















^0 








0/ 



(9) 



The adjoint of the mapping S is given by 



/Ao 



/AolA'-Y*(yo) 



V 









0\ 

Ail;f-Y*(yi) 

Yo 

Yi/ 



After a simplification of the primal and dual problems associated with (H, C, D), one obtains 
equivalent primal and dual problems as follows: 
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Dual problem 

minimize: l\\l^*{Y)\\^ + l\\wUy-') 
subject to: Y G Pd(Z). 



Strong duality 

To prove that strong duality holds for the semidefinite program above, it suffices to prove that 
the primal problem is feasible and the dual problem is strictly feasible. Primal feasibility is easily 
checked: one may verify that the operator 

/po \ 

pi 

Yo(|0o) 

Vo Yi(pi)/ 

is primal feasible for any choice of density operators po,pi G D{X). To verify that strict dual 
feasibility holds, one may consider the operator 

/Ao \ 

Ai 

I2 
Vo IzJ 

for any choice of real numbers Ao > || Yq(12) ||oo and Ai > || ^^(1^) ||oo. By Slater's theorem, 
strong duality follows. 

Optimal value 

For any fixed choice of po/Pi G D(A'), one has that the maximum value of the primal objective 
function 

lTr(X) + lTr(X*) 

subject to the constraint 

V X- T,{p,); 

is equal to F ( Yq (po ) a (pi ) ) / by the same analysis that was used to determine the primal optimum 
for the semidefinite program for the fidelity function. Maximizing over all choices of density 
operators po,pi G D(A:') gives Fmax(Yo/ ^i)/ which equals |||0|||i by Theorem|5l 



Primal problem 

maximize: ^ Tr(X) + ^ Tr(X*) 

/Yo(po) X \ 
subject to: > 

po,pi G D(A') 
X G L(Z). 
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3 A semidefinite program for the completely bounded trace norm from 
a mapping's Choi-Jamiolkowski representation 

In this section an alternate semidefinite program for the completely bounded trace norm is pre- 
sented. Whereas the semidefinite program from the previous section is obtained from a Stine- 
spring representation of a given mapping, the semidefinite program in this section is obtained 
from the Choi-Jamiolkowski representation of a given mapping. 

While the two semidefinite programming formulations are different, they are closely related. 
As for the semidefinite programs for the fidelity and the completely bounded trace norm in the 
previous section. Lemma |2] provides a key tool through which the semidefinite program given in 
this section may be analyzed. 

3.1 Choi-Jamiolkowski representations and the completely bounded trace norm 

Let A! = C" and y = C"', and assume that O G T( A", y) is a given mapping. The Choi-Jamioikowski 
representation of O is the operator /(O) S L(3^ (8) A!) defined as 

l<i,j<n 

An equivalent expression is 

/(<D) = ((D®lLw)(vec(lA')vec(l;r)*), 
where the vec-mapping is the linear mapping defined by the action 

vec(E,;y) = e,' ® ej, 

extended by linearity to arbitrary operators. 

One identity connecting the vec-mapping to the Choi-Jamiolkowski representation of a map- 
ping is the following one, which holds for all choices of A, B G L(A'): 

{\y ® A^) /(^) {\y®B) = {^® Il(a')) (vec(A) vec(B)*) . (10) 

Through this identity, an alternate expression for the completely bounded trace norm is obtained, 
as stated by the following theorem. 

Theorem 6. Let A" = C" and y = C", and let O S T(A', y) he a linear mapping. It holds that 
|||<I>|||i = max I 

Proof. By together with ^ it holds that 

= max{|| {\y ® A^)]{^){ty ®B) \\^ : A,B G L(A'), ||A||2 = ||B||2 = l}- 

By the polar decomposition, every operator X G L(A') with || X||2 = 1 may be written as X = ^/o^U 
for some choice oi a G D(A:') and II G By the unitary invariance of the trace norm, the 

theorem follows. □ 
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3.2 A semidefinite program from Theorem [6] 

The semidefinite program to be considered is specified by the triple (E, C, D), where 
is a Hermiticity-preserving mapping defined as 



V 



/Tr(Xo) \ 

Tr(Xi) 

Zo- 137^X0 

V Zi- 137^X1/ 



(11) 



and C e Herm(A' © ;f © (J^ ^ A") (J^ ^ A")) and D e Herm(C © C © (J^ O A") © (J^ O A")) are 
defined as 

/O \ 



and 




/(O) 
VO /(O)* J 

The adjoint of the mapping S is given by 

/Ao • • -X /AolA'-TryCyo) 
Ai ■ • _ 
• Yo • ~ 

V- • • V V 



D 



/I 












1 


















VO 





0) 




) 







°\ 


Trj; 


(Vl) 








) 




Yo 





) 










(12) 



After a simplification of the primal and dual problems associated with (E, C, D), one obtains 
equivalent primal and dual problems as follows: 



Primal problem 
maximize: ^(/(0),X) + i(/(0)*,X*) 



subject to: 



/I3; (g) po X ' 
1^ X* ty^pi, 

PQ,pi e D(A') 

X G L(3^^ A") 



> 



mmimize: 



subject to: 



Dual problem 

^l|Tr3^W|L + ^||Tr3;(yi)| 



/ Yo -/(0)\ 



> 



Strong duality 

Similar to the semidefinite programs discussed in the previous section, strong duality is easily 
established for the semidefinite program described above by the use of Slater's theorem. In fact, 
strict primal and strict dual feasibility hold for all choices of O; so that, in addition to strong 
duality, the primal and dual optima are achieved by feasible solutions in both cases. An example 
of a strictly feasible primal solution is 

/X 0\ 
0X00 
Z 

yo z/ 
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for 



X 



and 



dim(A') 

while an example of a strictly feasible dual solution is 

/A 0\ 
A 

y 
Vo y/ 



dim(A') ' 



for 



y 



+ 1 lyc^x and A = 1 



1/(^)1 



+ 1 dim(3;). 



Optimal value 

For any choice of density operators po,pi G D(<^), it holds that 



X* 



ty (g) pi 



> 



(13) 



if and only if 

X= {ly^^)K{ty®^) (14) 

for some choice of an operator K G L{y (S) X) satisfying 1 1 1 1 oo < 1, as follows from Lemma |2l The 
primal optimum is therefore given by 



sup ^{j{^),{ly(^^)K{ly(^^] 



sup 

Po,Pi 



ty®^i)]{^y{ty®^o] 



sup 

Po,Pi 



lloll 



1/ 



where supremums are taken over all K G L(A') with ||X||oo < 1 and pQ,pi G D{X), and where the 
last equality follows from Theorem [6l 



4 Remarks on the complexity of approximating optimal solutions to 
the semidefinite programs 

Suppose that (S, C, D) is an instance of one of the semidefinite programs described above, either 
for the maximum output fidelity characterization or the Choi-Jamiolkowski representation char- 
acterization of the completely bounded trace norm. It is natural to ask whether an approximation 
to the optimal value of this semidefinite program can be efficiently computed (under the assump- 
tion, let us say, that the complex numbers specifying S, C, and D have rational real and imaginary 
parts whose numerators and denominators are represented as integers in binary notation). 

From a practical viewpoint, algorithms employing interior point methods represent a sensible 
approach for computing the optimum value of these semidefinite programs PA li95l [dK02L The 
CVX software package [GB09J for the MATLAB numerical computing environment allows one to 
solve these semidefinite programs efficiently with minimal coding requirements. 
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For the sake of obtaining rigorous statements about the polynomial-time solvability of the 
semidefinite programs (and perhaps not much more than that), the ellipsoid method is a more at- 
tractive alternative, applied specifically to the dual formulations of the semidefinite programs. 
When considering the applicability of the ellipsoid method, it is helpful to consider the following 
set, for V C Herm(3^) denoting the dual feasible set of (H, C, D) and e > being a positive real 
number: 

V° = {y e Herm(3^) : Y + H G P for all H G Herm(3^) satisfying ||H||2 < e}. 

Intuitively speaking, contains every operator in the interior of the dual feasible set that is not 
too close to the boundary of that set. 

It has already been demonstrated that 1)° is nonempty for some choice of e for each of the 
semidefinite programs, in the discussions of strong duality in the two previous sections. To argue 
that accurate approximate solutions to the semidefinite programs can be obtained by the ellipsoid 
method, a sufficiently large lower bounds on the value of £ for which 1)° is nonempty is needed. 

For the semidefinite program for the maximum output fidelity characterization of the com- 
pletely bounded trace norm, presented in Section|2l the adjoint of the mapping H is given by 
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for 

Ao = ^ + ||YS(l)||oo and Ai = l + ||Y5[(1) 
is a specific example of a strictly dual feasible solution satisfying 
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A calculation reveals that for H G Herm(C (B C ® Z ® Z) satisfying 

||H||2 < ^min{||YS||oo^ l} 

it holds that ||H*(H)||oo < 1/2. As Y* and Y* are positive, it holds that ||Y*||oo = ||Y*(l)||oo and 
II Y^ ||oo = II Yj (1) ||oo, from which it follows that V" is nonempty for 

1 

^~ 4(l + ||Y*(l)||oo+||Y*(l)||oo)' 
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For the semidefinite program for the completely bounded trace norm presented in Section |3l 
based on the Choi-Jamiotkowski representation, the adjoint of the mapping H is given by 
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is an example of a strictly dual feasible solution satisfying 
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{y(SiX)) satisfying 
1 



2 < 



it holds that ||E*(H)| 



2dim(J^) 

< 1, from which it follows that 1?° is nonempty for 

1 



2dim(3^) 

In both cases, the lower bound on the value of e for which is nonempty is polynomial in 
the input data and efficiently computable. 

One also requires an upper bound on the size of an optimal, or near optimal, dual feasible 
solution. For the semidefinite program based on the maximum output fidelity characterization of 
the completely bounded trace norm, every dual feasible solution is positive semidefinite, and for 
approximate solutions it is sufficient to consider only those dual feasible solutions whose trace is 
at most 

R = Hoc + 11^^(1) ||eo+2dim(Z). 

For the semidefinite program for the Choi-Jamiotkowski representation characterization of the 
completely bounded trace norm, every dual feasible solution is again positive semidefinite, and 
an optimal solution cannot have trace larger than 

R = 2||/(0)||eodim(;f)dim(3;). 

The trace of every positive semidefinite operator serves as an upper bound on that operator's 
Frobenius norm, which implies that the above quantities also upper-bound the Frobenius norm of 
the set of dual feasible solutions that are worthy of consideration. 

As is described in detail in IIGLS93II for a significantly more general setting, and summarized in 
|Lov03B for the semidefinite programming setting, the bounds e and R above allow one to conclude 
that an algorithm running in time polynomial in the input size and log(l/^) can approximate the 
optimal value of the semidefinite programs discussed above to within accuracy S. 
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